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In this paper we investigate the evolution of binary neutron stars, namely, their magnetic field, 
spin, and orbital evolution. The core of a neutron star is considered to be a superfluid, supercon- 
ductor type II. Flux expulsion of the magnetic field out of the core of a single neutron star has been 
discussed by previous authors. However, the evolution of the core magnetic field is substantially 
different for a binary neutron star. While for a single neutron star the fiuxoids of the proton su- 
perconductor always move outward through the core, in a binary neutron star in accretion phase 
fluxoids move back into the core. The subsequent increase of the core magnetic field results in 
the increase of the surface magnetic field. We ask whether within the framework of this model 
the formation of millisecond pulsars (MSPs) is possible. We show that despite the increase of the 
core magnetic field, MSPs are formed in this model. The evolution of neutron stars with various 
orbital periods, magnetic fields, spin periods, and other parameters are numerically investigated. 
The equation of state of the neutron star, initial orbital period of the binary, and other parameters 
of the binary have substantial effects on the evolution of period vs. magnetic field. 
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I. INTRODUCTION 



The evolution of the magnetic field of neutron stars has 
been widely discussed in many papers [H-Q- The present 
study is mainly addressed to the coupled spin, magnetic 
field, and orbital evolution of a binary neutron star. In 
this paper we investigate motion of fluxoids into and out 
of the superfluid, superconductive core of a neutron star 
and the subsequent evolution of the magnetic field in the 
crust. To derive the magnetic out-flux and in-flux rates, 
we consider various forces which act on fluxoids in the 
interior of a neutron star ^ and 0. The evolution of a 
binary star is affected by accretion from the stellar wind 
matter of its companion by various effects. These effects 
were discussed in ^5] under the assumption that magnetic 
field is confined within the solid crust of a neutron star. 
In this paper we extend the results of [2| (for isolated 
neutron stars) to binary neutron stars, taking into ac- 
count the various effects of the accretion discussed in Q 
and We show that accretion significantly modifies 
the evolution of the magnetic field. Namely, within the 
framework of this model, the magnetic field of a binary 
neutron star can increase, while magnetic field of a single 
neutron star always decreases. 

As a neutron star evolves, the strength of its surface 
magnetic field and spin period change. The Ps — Bs dia- 
gram is usually used to present distribution of pulsars of 
different surface magnetic fields Bs and spin periods Pg 
(here 's' stands for spin, not surface), or their evolution. 
Observed pulsars are distributed within two separate re- 
gions in Ps — Bs diagram. MSPs are fast spinning stars 
{Ps < 10—20 ms). Their surface magnetic fields Bs is less 
than 10^" G. Other pulsars (normal pulsars) have higher 
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magnetic fields and longer spin periods. Almost 50% of 
MSPs arc in binaries, while only 2% of normal pulsars are 
found in binaries. It is believed that newly born neutron 
stars have high surface magnetic fields and short spin pe- 
riods. The spin period of single neutron stars increases 
during time, and their magnetic field decays. It has been 
suggested that MSPs are old binary neutron stars. Ac- 
cording to the standard picture Q, @, binary neutron 
stars pass through four evolutionary phases: the ejector, 
the propeller, the accretion from the wind of the com- 
panion, and the accretion resulting from the Roche-Lobe 
overflow. This four-phase model explains the low mag- 
netic field and fast spin period of the millisecond pulsars, 
the so-called recycled binary neutron stars. 

The mechanism of generation of the magnetic field 
of neutron stars is not totally understood. However, it 
has been suggested that magnetic field might penetrate 
through both the core and the crust of a neutron star 
[loj . However, as we have mentioned in our model, the 
evolution of the core magnetic field and subsequently the 
surface magnetic field are significantly different for binary 
and single neutron stars. In a binary neutron star, the 
magnetic field increases during the accretion phase due to 
the in-fiux of the magnetic field into the core. This raises 
the intriguing question whether MSPs can still be formed 
within the framework of such models. We shall illustrate 
that if one includes the core magnetic field, there still ex- 
ists the possibility for the formation of MSPs in binaries. 
However, not all initial values of the core magnetic field 
lead to the formation of MSPs. 

As it was discussed in Q, the orbital separation of a 
binary affects the rate of change its neutron star spin pe- 
riod. However, the orbital separation also evolves during 
the evolution. We present the importance of the consid- 
eration of the orbital evolution on the final stage of a 
neutron star. 

The decay rate of the magnetic field in the crust de- 
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pends on the crustal conductivity. The conductivity de- 
pends on the equation of state (in the crust) and on the 
temperature, and hence, on the age of the neutron star. 
Therefore, the equation of state in the crust affects the 
evolution of a neutron star in the Ps — Bg diagram. Ac- 
cording to 5] , the evolution of neutron stars with a stiff 
equation of state like the Friedman-Pardharipande are in 
agreement with the Ps — Bg distribution of the radio pul- 
sars. In a model for the evolution of the magnetic field, 
where the magnetic field is not confined to the crust such 
as the model of Konenkov and Geppert , the relation 
between the equation of state and the magnetic filed de- 
cay is apparent as well. Here we did not consider the 
effect of various equations of the states. However, we 
propose that a theoretical model such as this model can 
be compared with the observational Ps — Bs distribution 
of pulsars to rule out some of equations of the states. In 
fact, the numerical result of this paper can be extended 
to various equations of the states and cooling scenarios. 
The population synthesis is a necessary calculation for 
setting a limit on possible equations of states within the 
framework of the present model. These calculations can 
be done in future works. However, the main goal of this 
paper is to remind a reader of the importance of the ex- 
tension of the results of 2] for single neutron stars to 
binary neutron stars and to present the result of such an 
extension in a few examples. In this paper we investigate 
the history of recycled binary neutron stars taking into 
account the various effects associated with the supercon- 
ductive, superfluid core. We show that the presence of 
these effects modifies the evolution of the surface mag- 
netic field substantially. The surface magnetic field of 
a binary neutron star with superfluid, superconductive 
type II core may increase during the accretion phase. 
We ask whether this alters the formation of MSPs. We 
show that despite the increase of the core magnetic field, 
MSPs are formed. In Sec. II we present a review of the 
standard evolution of binary neutron stars. The model 
for the evolution of the surface and core magnetic field is 
explained in Sec. III. In Sec. IV we present our model 
for the evolution of binary neutron stars combining the 
general features of the Sec. II and Sec. III. In Sec V we 
give an overview of the corresponding numerical results. 
Throughout this paper we discuss the evolution of the 
orbital separation of a binary a(t), the surface magnetic 
field of a neutron star Bs{t), the core magnetic field of 
a neutron star Bc{t), the mass of a neutron star Mn{t), 
and the mass of the companion M2(t). Therefore, any 
quantity depending on the above quantities also depends 
on time implicitly. Throughout this paper we avoid re- 
peating the dependence on time of various quantities. 
However, the reader should be aware of this dependence 
unless otherwise stated. 

The main features of the model such as: 1) the in- 
crease of the strength of the magnetic field in the accre- 
tion phase, 2) the independence of the final strength of 
the surface magnetic field from the initial values of the 
core or surface magnetic fields, 3) the crucial role of the 



TABLE I: Values of mass M2, lifetime T2, mass- loss rate M2 
and wind velocity V-w for low-mass (LM), intermidiate-mass 
(IM), and high- mass (HM), and Be- type companion stars 
(CS). 



cs 


M2 


T2 


- log M2 


v„ 




(Mq) 


(yr) 


(Mgyr-i) 


(kms-^) 


LM 


0.8 - 1.0 


3 X lO'' - 10"' 


12 - 16 


400 - 800 


IM 


2.5-7 


5 X lO'^ - 7 X 10** 


9-11 


500 


HM 


> 9 


1.2 X 10^ - 3 X 10^ 


6-10 


600 - 700 


Be 


9 


3 X 10^ 


8-10 


200 



minimal spin period of the pulsar on the magnetic field 
evolution, are independent of the chosen parameters or 
the equation of state. 



II. EVOLUTION OF NEUTRON STARS IN 
BINARIES 

The weak magnetic field and fast rotational spin of 
MSPs suggest the idea of enhanced magnetic field decay 
in accreting neutron stars. According to the standard 
picture, MSPs are old neutron stars spun up by accre- 
tion of matter from their companion in binary systems. 
The accretion also heats up the crust of a neutron star. 
Hence, accretion enhances the decay of the magnetic field 
due to the decrease of the crustal conductivity. There- 
fore, in order to track the formation of MSPs in binaries, 
we consider a model for orbital, spin, and magnetic field 
evolution of a newly born neutron star of mass M„ and 
radius i?„ in a binary with initial orbital period Poi- The 
companion has mass M2 and radius i?2. In the study of 
binaries one might consider different companions. Table 
I presents the values of different parameters that can be 
adopted for different companions. In the case of a low 
mass companion, wind accretion does not produce X-ray 
emission. However, the long main sequence lifetime of the 
companion allows for a long evolution in different phases. 
The companion loses mass in the form of a spherical uni- 
form stellar wind at rate M2. The orbital separation of 
the binary and consequently orbital period, change dur- 
ing the evolution due to the accretion, mass loss from the 
binary, and two other sinks of orbital angular momentum, 
namely, magnetic braking and gravitational radiation. In 
the case of a low mass companion, the orbital separation 
of the binary remains almost constant during the stellar 
wind evolution. However, for short initial orbital periods 
orbital evolution can alter the evolution of the surface 
magnetic field dramatically at the end stages of the evo- 
lution. Fig [in] shows that including the effect of orbital 
evolution can change the final value of the final surface 
magnetic field Bsf up to one order of magnitude. There- 
fore, we do take into account the effect of the loss of the 
angular momentum from the system. 
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In a spin- induced field expulsion scenario (SIF) [Tl|, 
the evolution of the spin period and the magnetic field of 
a neutron star would be coupled. Our model is different 
from the original SIF model, as it will be explained later, 
but the coupling feature still persist in this model. 



A. The binary neutron star vs. phases 

The standard evolution of neutron stars in close bina- 
ries is described in four successive phases: the ejector, 
the propeller, the wind accretion, and the Roche Lobe 
overflow [H']. In principle these phases might not be con- 
sidered successively. A neutron star might return to a 
previous phase. However, the conditions for entering a 
phase can be different from that of returning back to a 
phase. This also brings a lot of complications on the tem- 
perature of the neutron star versus time. For example, 
for the transitions between the propeller and the accre- 
tion phase, a fight between the free cooling of the neutron 
star and heating due to the accretion arises. These con- 
ditions have to be well studied before being included in 
any numerical computations. One should be aware that 
including such possibilities will change our picture of the 
evolution drastically. For example, a neutron star might 
reach almost a constant value of magnetic field and spin 
period (before reaching an appropriate value for a MSP) 
while oscillating between accretion and propeller phase. 
This might alter the formation of the MSPs. Here, We 
consider that these phases occur successively. Therefore, 
all our results and conclusions depend on this assump- 
tion in the first place. Here, we give a brief overview of 
these phases. 

Ejector: A neutron star strong radiation pressure keeps 
the companion stellar wind matter away from its magne- 
tosphere. Hence, a neutron star does not experience any 
influence from its companion and evolves like an isolated 
neutron star (single neutron star) . During this phase the 
magnetic field of a neutron star decreases and its spin 
period increases. We assume that during this phase a 
neutron star cools down according to the standard cool- 
ing scenario [l^l • 

Propeller: Due to the decrease of magnetic field the 
pulsar radiation pressure decreases. Consequently, at 
some stage the neutron star radiation pressure is not suf- 
ficient enough to keep the companion stellar wind matter 
away from its magnetosphere, and the stellar wind pen- 
etrates into the neutron star magnetosphere. The mag- 
netosphere of the neutron star ejects the infalling wind 
matter transferring some part of the neutron star angu- 
lar momentum to the wind. Therefore, the neutron star 
spins down. The accretion torque produced by the cap- 
ture of the plasma from the stellar wind depends on the 
detailed properties of the outflow. A steady Keplerian 
disk may form from the accretion flow outside the mag- 
netosphere. Spherically symmetric radial infall may also 
occur for these systems [l3-[3- We use a coefficient C 
to count for the details of the flow In a more ac- 



curate model, numerical simulations of the flow should 
be included as well. In this phase thermal and magnetic 
field evolution of neutron stars do not differ from that of 
single neutron stars. 

Wind accretion: Finally, the magnetic field of a neu- 
tron star decreases sufficiently enough that accreted mat- 
ter falls onto its surface. In this phase the thermal, mag- 
netic field, and spin evolution of a neutron star differs 
completely from that of a single neutron star. Accretion 
with an accretion rate M2 = 10~^° — lO~^M0yr~-'^ or 
higher heats up the neutron star to an accretion temper- 
ature Tac = (2 - 3) X lO^K and 0. 

The Roche lobe overflow: This phase begins when the 
companion star fills its Roche lobe, so the mass trans- 
fer is strongly enhanced. Accretion resulting from Roche 
lobe can last as long as (1 — 5) x lO^yrs in low mass 
binaries. We assume that a neutron star enters Roche 
lobe overflow when cither the companion leaves its main 
sequence (for the life time of different mass stars refer 
to Table I) or the radius of the companion R2 equals to 
the Roche lobe radius Rl- This is possible to happen 
because we are taking into account the loss of orbital 
angular momentum of the binary and decrease in the bi- 
nary radius a; however this is very unlikely for a low main 
sequence companion . A measure of the Roche lobe 
radius is given by the radius Rl oi a sphere which has 
the same volume as the lobe (the critical equi-potential 
surface which passes through the inner Lagrangian point 
Li). A simple expression for the Roche lobe radius Rr 
is given only in term of the ratio of the mass of the two 
stars q = Ah/Mn [Hi, 

Rr = a(0.38 + 0.21ogg) (0.5 < g < 20) 



R 



R 



a(0.462 



) (0<g<0.5), (1) 



where a is the orbital radius of a binary (|16p . 

We shall perform our numerical calculations only for a 
low main sequence companion. The computations are 
performed for a lifetime equal to the expected main- 
sequence lifetime of the companion star. Computations 
are stopped before the Roche-lobe overflow. Note that 
there is no difficulty in extending the results for other 
companions and beyond Roche-lobe. 

To define the phase of a binary neutron star, we com- 
pare the accretion radius i?a, the radius of magneto- 
sphere and the corotation radius Rco- Let us first 
define each of these radii. 

Stellar wind matter of impact parameter less than 
the so-called accretion radius, is absorbed by a neutron 
star gravitational field. The accretion radius is given by 



where 



Ra 



2G'M„ 



(2) 



(3) 
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is the orbital velocity of a neutron star and Vw is the ve- 
locity of the wind from its companion [19l] . The wind ve- 
locity changes when the distance from the star increases. 
The wind is believed to be accelerated to a typical ter- 
minal velocity of a few times the escape velocity at the 
stellar surface [131 . However, we consider a constant 
wind velocity assuming that the wind velocity close to 
a neutron star has already achieved its terminal velocity. 
Therefore, is considered as a constant parameter in 
our model. 

Not all the mass lost from companion in the form of 
stellar wind is accreted by a neutron star. The accretion 
rate reads 



(4) 



where C'oo is the speed of sound in the stellar wind mat- 
ter, Vrei = VmI + K) I ^■iid r is the distance from the 
companion [13]. We assume that M2 is constant during 
the life-time period of the companion. The energy radi- 
ated by a rotating neutron star in the form of relativistic 
particles or electromagnetic radiation originates at the 
characteristic distance of the so-called light cylinder ra- 
dius 



(5) 



where c is the speed of light and fig — 2tt/Ps is the 
spin angular velocity of a neutron star. The stellar wind 
matter is a highly conductive plasma, and at distances 
r < Ri the electromagnetic field is static. Therefore, 
the wind can not penetrate through the field, but it can 
squeeze the field. The radius of the magnetosphere Rm is 
defined where the pressure of the magnetic field around 
a neutron star and the pressure of the matter balance 
each other. At a radius r > Ra the pressure of the stellar 
wind (r, t) , is given by the following equation 



p — 



M2K, 



47r(a - 



(6) 



At a radius r < Ra the stellar wind pressure grows as 
r~2 J and it is given by 



Py_, =r-2Ma 



87r2 



(7) 



If at r = Ra the pressure of the dipolar magnetic field of 
a neutron star is smaller than the pressure of the stellar 
wind matter, the radius of magnetosphere reads 



[BsR„ 



(8) 



where Bg is the strength of the magnetic field on the 
surface of the neutron star (3] . For the case of a disk ac- 
cretion, R„i is similar. However, both smaller and larger 
values have been given by different authors |37(| and (U] . 



If the radius of magnetosphere is larger than Ra, it can 
be found from the equation 



2\M2 



rel 



Rl, + (i?™ - a) = 0. 



(9) 



The magnetosphere of a neutron star is bounded by the 
light cylinder of radius Ri Q. 

To find the condition for the transition of a pulsar from 
the ejector to the propeller phase we first introduce the 
radiation pressure Prad, pressure due to particles and ra- 
diation field emitted by a neutron star at the distance r 
from the star, 



Prad — 



rad 



(10) 



where Lrad is the total rate of the energy loss by the 



neutron star 



2 
3^ 



RtB^Clt. 



(11) 



Rrad is defined as the radius where the radiation pressure 
and the pressure of the stellar wind matter balance each 
other. Here we have two cases: 

(1) Ri < Rrad {Prad > Pgas at r = Ri): & ueutrou 
star will be in the ejector phase or the (most probably 
obscured) active radio pulsar phase. 

(2) Rrad > Ri {Prad < Pgas at r = i?; ) : the accretcd 
matter makes contact with the magnetosphere and may 
cause the spin-up or spin-down of a neutron star. When 
the matter reaches the magnetosphere it interacts with 
the neutron star; hence the angular momentum of the 
neutron star changes by the following equation 



Lg 



-(MaRMr^R^), 



(12) 



where C is an efficiency factor included to take into ac- 
count the details of the geometry of the interaction (C = 1 
corresponds to the case of disk accretion), and Vd(i?m) = 
Vco{Rm) — Vkep{Rm) is the difference between the coro- 
tation velocity of stellar wind matter, Vco{Rm) = ^sRm, 
and Keplerian speed of the stellar wind matter around 

the neutron star, Vkep{Rm) — \J , both evaluated at 
distance i?„i from the neutron star. 

Define the radius of the corotation surface in the fol- 
lowing form 



Rc 



(13) 



When Rco < Rm, a neutron star is in propeller phase 
and the spin period of the neutron star decreases, but 
when Rco > Rm the star is in wind accretion or spin-up 
phase. The amount of spin angular momentum added or 
extracted from the system is calculated assuming that: 
1) In the ejector phase all the stellar wind of the com- 
panion leaves the system. 2) In the propeller phase all 
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the accreted matter is expelled by the neutron star at its 
magnetosphere with an angular velocity equal to that of 
the neutron star. 3) In accretion phase all the accreted 
matter enters the magnetosphere. In all these cases the 
fraction of the stellar wind that does not accrete to the 
neutron star leaves the binary with the angular momen- 
tum equal to that of the companion. 

The spin-up of a neutron star due to the accretion is 
limited to an equilibrium period given by 



eg 



2ab: 



(14) 



where Bs9 is the surface magnetic field in the unit of 
lO^G, RnQ is the radius of the neuron star in the unit of 
lO^cm, MEdd is the maximum possible Eddington limit 
for the accretion rate given by 



MEdd 



1.5 X 10"^i?„6 



M, 



(15) 



The spin-up line in the pulsar Bs — Ps diagram is given 
by Peq{MEdd), where Ma is set equal to MEdd- In our 
calculation we set this value as the minimum value of the 
period that a neutron star can get to due to the accretion. 
This limit plays a crucial role in our model. This limit 
controls the increase of the core and the surface magnetic 
field of a neutron star in the accretion phase. 



B. Orbital evolution 

We have a binary system composed of a neutron star 
of mass Mn and its companion of mass Af2- We assume 
that both stars are in circular orbit around the center 
of mass. Then orbital angular momentum of the neu- 
tron star and its companion is given by L„ = MnRj^flo 
and L2 = M2R2^o, respectively. Here, ilo is the orbital 
angular velocity of each star around the center of mass, 



GM 



(16) 



where a is the orbital separation (orbital radius) of the 
binary and M = Af„ + A/2 is the total mass of the binary. 
For our binary system the orbital angular momentum , 
Lo = Ln + L2, is given by 



, , , , / Ga 
Lo = MnM2\ — . 

V M 



(17) 



Now, we assume that the companion is losing mass in 
the form of stellar wind. Some fraction of this mass will 
leave the system. Also, there are some other sources 
of the angular momentum loss from the system such as 
gravitational waves or magnetic braking. All these effects 
cause the dependence of a, M, M„, M2, and Lq to time. 
The change in a, affects the rate of change in the spin 



period of a neutron star Pg . From ((T7)) the orbital radius 
evolves according to the equation 



a 
a 



where 



Here, 



Lo 



M 
M 



L 



,Mn 



,M2 
'M2'' 



MB 



Law- 



(18) 



(19) 



Lmb = -0.5 X 10-"'f-'k'M2 Rl K g cm^s-^ (20) 

is due to magnetic braking (23| , f2o is the orbital angular 
velocity of the binary, / = 1.78, and k = 0.1 J^. The 
rate of the loss of the orbital angular momentum from the 
binary due to the gravitational radiation Lqb is given by 



Lgw = —32- 
Finally we have 



^G^MlM^M^ 
5c^a5 



(21) 



(22) 



where Lesc is the rate of the loss of the angular momen- 
tum from the binary due only to the escape of stellar 
wind matter \7] and /? is the fraction of the angular mo- 
mentum of the companion star that escapes from the 
binary when the rate of total mass loss from the binary 
is M = Ma — M2- In the accretion phase /3 is given by 
(see [2^ for details of the derivation) 



(3 = 1 



a — 1 



RmMi\Vd\ 
G-2a-2{t)M^ 



(23) 



and in the propeller phase by 
M2 

Mn 




Rn 



M 
M^ 



(24) 



Here, a = (1 — Ma/M) is the fraction of the stellar wind 
matter which leaves the binary without any interaction 
with the neutron star. The rest interacts with the magne- 
tosphere and changes the rotational angular momentum 
of the neutron star. In terms of the above parameters, 
(|17p can be written in the following form: 



,A/2 aM2 Mn 
"^M:" 2M -"^Af 



M2 
M2 



L 



MB 



Lo 



(25) 



We assume that the initial values of the spin period 
and of the magnetic field of a neutron star born in a 
binary are independent of its orbital period. 
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III. MAGNETIC FIELD EVOLUTION 

We assume that initially the magnetic field of a neutron 
star penetrates both into its core and crust. We do not 
discuss the mechanism of the creation of this magnetic 
field. We assume that in the crust the field is maintained 
by the currents. We restrict ourselves to a consideration 
of the decay of a field which is assumed to be dipolar 
outside the neutron star. The evolution of the magnetic 
field B in the crust of the neutron star is given by [1^ 
and 1221 



9B(r,t) 
dt 

where 



-V X (^^V X B(r, t)j + V X (v X B(r, t)) , 

(26) 



Ma 

V = VBr = ; —Gr, 



(27) 



is the velocity of the material movement, and a is the 
electrical conductivity in the crust. The velocity of the 
material movement is zero in the ejector and propeller 
phases. The magnetic field B is obtained from a vector 
potential B = V x A. Due to axial symmetry, in spher- 
ical coordinates, A = {0,0, A^{r, 9,(1))). Introducing 



we have 



S{r,t) sinO BsoRnS{r,t) sint 



(28) 



Brir,t) 
Bg{r,t) 



2 sine ds{r,t) 



(29) 



where Bgi is equal to the value of the surface magnetic 
field at 9 = 7r/2, where the value of the surface mag- 
netic field is maximum. |B|(r = Rc) = -BsiVS cos^^ -I- 1. 
Therefore, from equation (I26p we get 



dS{r,t) 
dt 



Ana 



d^Sjr, t) 2S{r, t) 

Qj.2 j,2 



+ v{r, t) 



dSjr, t) 
dr 



(30) 

Boundary condition on the surface of a neutron star reads 



dS{r, t) 



dr 



S{Rn,t)=0. 



(31) 



The above boundary condition ensures that the magnetic 
field is dipolar outside a neutron star. 

Conductivity comes mostly from the scattering of the 
electrons on phonons and on impurities [29i|. The con- 
ductivity due to the scattering on phonons apt is found 
from 



1.21 X 10 



28 • 



J12 



rp2 , 



0.084r2 



(32) 



where T is the temprature which changes with time ac- 
cording to the standard cooling scenario (except in the 
accretion phase) , and To is the Debye temperature 



(33) 



where Z and A are the atomic number and the mass 
number of the dominant element in the crust of a neutron 
star. For a table of the dominate elements at different 
densities in the BPS equation of state refer to Here, 



Pf 



vrieC 



(34) 



where Pp is the Fermi energy, m,, is electron mass, c 
is the speed of light, x can be written in terms of the 
density in the following form 



h 



ArrLn 



(35) 



where me is the mass of electron, rup is the mass of 
proton, and h is reduced Plank constant. The phonon 
conductivity decreases when T increases, therefore addi- 
tional heating due to the accretion accelerates the mag- 
netic field decay. At very low temperatures aph is very 
large. The conductivity due to impurities reads 



imp 



= 8.53 X 10^ 



•^^^imp Q 



(36) 



where Ai„ip is the Coulomb logarithm for impurity scat- 
tering (at p > lO^gcm"'^ one has A^^p = 2), and Q 
is the parameter characterizing the concentration and 
charge of impurities, known as the impurity parameter 
[2^. Here we assume that 0.001 ^ Q ^ 0.1 and that Q 
is constant throughout the crust for the whole period of 
the evolution. We assume that at very low temperatures 
T < lO^K, the conductivity is given by only the scatter- 
ing on the phonons. A binary neutron star in most cases 
will not reach to such temperatures due to accretion (ac- 
cretion starts before the cooling makes the neutron star 
so cold). 

Finally, the total conductivity in the crust, a, is given by 



1 
a 



1 

<7ph 



1 



(37) 



Magnetic field evolution of the neutron star is sensitive 
to the temperature due to the dependence of the conduc- 
tivity on the temperature. Accretion effects the crustal 
conductivity substantially. We assume that temperature 
is uniform throughout the crust. For more exploration of 
the current model one has to consider the dependence of 
the temperature to the density. 

We assume that in the core of a neutron star protons 
form a superconductor of type II, and the magnetic fiux 
is concentrated in the array of magnetic fiux tubes called 
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fluxoids. The evolution of the magnetic field in the su- 
perconductive core of a neutron star is given by 



dB^r, t) 

dt 



V X (vp(r,<) X Bc(r,t)), 



(38) 



where Vp(r, i) — Vp{r,t)er is the radial velocity of flux- 
oids. 

The core of a neutron star consists of neutrons in a su- 
perfluid state and protons in a superconductor state. In 
steady state vortices in the core of a neutron star coro- 
tate with the superfluid bulk matter at the rate fi. When 
a star spins down at a rate f2 , the superfluid bulk matter 
rotates at a faster rate, fib. We assume that if fi is kept 
constant there is a constant difference (lag) between the 
angular velocity of the superfluid bulk matter and the 
neutron vortices = fif, — fis > 0. The radial velocity 
of the vortices at the core-crust boundary is determined 
from the conservation of angular momentum [l| and (llj 

where is the radius of the core. When a neutron star 
spins down Ps > 0, and vortices move outward, so one 
has v„ > 0. When the neutron star spins up, P, < and 
vortices move inward to the core, so one has v„ < 0. For 
a single neutron star Pg is always positive and vortices 
always move outward from the core. In a binary neu- 
tron star in the accretion phase, the neutron star spins 
up. Therefore, vortices start moving back into the core. 
In the original (SIF) model jllj the pinning between the 
neutron vortices and the proton fluxoids was suggested 
to result in the expulsion of the magnetic flux out of the 
core, into the crust of a neutron star at a rate equal to 
the spin-down rate of the star. However, due to the ex- 
isting resistance against the motion of fluxoids and vor- 
tices, fluxoids and vortices move with different velocities. 
There are several forces acting on fluxoids driving them 
outward to the conductive crust of a neutron star, where 
we have an ohmic decay (pO|. Therefore, to find the rate 
of the out-flow or in-flow of the magnetic flux we consider 
various forces which act on the fluxoids in the interior of 
a neutron star such as the force due to their pinning in- 
teraction with the moving neutron vortices, the viscous 
drag force due to the magnetic scattering of electrons, 
and the bouyancy force. The resulting motion of flux- 
oids in a binary neutron star is not trivial. However, 
there exist a possibility for the outflow of the fluxoids. 
Bc{Rc, t) gives us the appropriate inner boundary condi- 
tion (at the bottom of the crust) for equation ([50)1 . The 
pinning per unit length exerted on a fluxoid is given by 

m 



f„(i?c, t) 



2<^opR^ns{t)uji{t) ^ 

Bc{Rc, t) 



dyn, 



(40) 



where $o = 2 x 10 cm^ is the quantum of the mag- 
netic flux carried by the fluxoid, p is the density in the 



core-crust boundary of a neutron star, and Bc{Rc,t) is 
the strength of the magnetic field in the core-crust bound- 
ary in units of G Pinning can be negative or positive. 
The maximum toi which can be sustained by the pinning 
force defines the maximum force which can be exerted on 
fluxoids by vortices. This maximum value of uji is known 

as UJcr ISI] 



1.59 X 10" 



^B^g^ rad s 



(41) 



where Bcs is the strength of the core magnetic field in 
the units of lO^G. 

The bouyancy force per unit length exerted on fluxiod 
reads 



)Ouy; 

m 



Note that 



Xp = 1.315 X 10"^^ 



cm, 



(43) 



Xp is the London penetration depth, where ^ > V2 for 
a superconductor type II. The bouyancy force is always 
positive, that is it is always directed outward. The drag 
force per unit length exerted on fluxoid is given by 



URc,t) 



3tt n. 



e^<i>lvp{R„t) 



64 



EpXp 



dyn. 



(44) 



where is the number density of electrons, Ep is the 
electron Fermi energy. Here we use the values rie = 3 x 



10' 



,36 



1.7 X lO'^** cm-'^, and Bp 



Mev m 



When Vp > 0, that is when fluxoids are moving outward 
f„ < 0. 

It was shown in [sij that the migration of fluxoids de- 
pends on the rate of the URCA process due to the fact 
that it is dependent on the beta equilibrium and the rate 
it is restored. Here we consider the model proposed by 
[3^ . since we are closely following the formalism devel- 
oped by and We are interested in comparing our 
results with 0. Namely, we want to compare the sin- 
gle and binary stars while we keep the formalism of the 
magnetic field evolution the same. One may improve the 
model considering the fluxoids drag derived in [31] rather 
than dH. 

Assuming that the magnetic flux in the core is uniform, 
we can find the total forces exerted on fluxoids [2,] 



4i?. 



(45) 



where 4/3i?c is the mean length of the fluxoids, and Np = 
ttR^Bc/^q is the number of fluxoids. The net power due 
to the forces is equal to the Poynting flux through the 
surface of a neutron star core 



E 

fluxoids 



/ {fb + fn + fv)Vp dl = --^ / [Ex B].dacore, 

J J Score 

(46) 
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where the integration on the right hand side is done over 
the surface of the core, the normal vector of this surface 
is being directed inward [2]. This equation can also be 
written in the form [3| 



F„+F6 + F„(vp)+Fe = 0. 



(47) 



The above equation gives us the radial velocity of the 
fluxoids, Vp. From equation (j46p one can find the appro- 
priate form of Fcrust- We have 



„ ^ 47r , 
V X B = —J, 

c 



(48) 



therefore, we find 



E 

and using ([26 

Fc{Rc,t) 



c I sm 



r 



d^S ^ 25 



(49) 




(50) 



where ctc is the conductivity at the core-crust boundary 
and Vc = Ma/ {4:TtRIpc)- The second term is 10^" times 
smaller than the first term, and is zero except in the 
accretion phase. We neglect this term. Making the same 
assumptions as we have made deriving equation (j30p . we 
can write equation (|38p in the from 



dS{r, t) 



-Vp{r, t) 



dS{r,t) 



(51) 



dt ' ' dr 

We assume that the core magnetic field is homogeneous 
therefore 

Bcity 



Equations ([5 



Sir,t) 

and ([5T|) give 

Vp{r,t) = j{t)r, 



(52) 



(53) 



Making the assumptions above, we have 



5(r = i?„t) = 5(r = i?„io)exp(-2 / 7(0^^0, (54) 

Jo 

where 5(i?c,io) is the value of the stream function in 
the core-crust boundary at the initial time to- When 
fuxoids are transported out of the superconductive core 
(vp > 0), the value of the stream function S at the core- 
crust boundary decreases. When Vp < 0, fluxoids move 
back into the core and the magnetic field in the core in- 
creases. (HH) gives us the inner boundary condition 0] for 
equation ([50]) . S{Rc,t). The increase of S{Rc,t) (when 
Vp > 0), leads to the increase of the surface magnetic 
field of a neutron star. From ([52l and |44)) Fc is given by 
(the second term is neglected) 



IV. DESCRIPTION OF THE MODEL 

We assume that t > 10* yrs of the evolution of binary 
neutron stars will result in the formation of MSPs. Ac- 
cretion of matter from the stellar wind of the companion 
affects the magnetic field evolution due to the change of 
the crustal conductivity, the orbital evolution, and most 
importantly the spin evolution of a neutron star. Here, 
we consider the coupled spin, mass, magnetic field, or- 
bital separation, and orbital period evolution of a neu- 
tron star in a binary. Namely we combine the magnetic 
field model discussed in Sec. Ill to the standard spin, 
and orbital evolution of a neutron star. We also take 
into account the cooling of the neutron star. Note that 
the magnetic field model we consider here was previously 
considered only for single neutron stars. 

In the ejector, propeller and accretion phase we inte- 
grate a set of the differential equations (to be explained 
later), to find the evolution of quantities M„(t), a{t), 
Ps{t), s{t), Bs{t) and Bc{t). To integrate equation ((50)) 
we have to know the initial stream function and the in- 
ner (j44p and the outer boundary conditions. 7 in (j54p is 
proportional to the velocity of the fluxoids. The veloc- 
ity of the fluxoids, Vp, is given by equation ([17)) . How- 
ever, in this equation w; is unknown as well (since Fn 
depends on wi). The following strategy helps us on the 
derivation of Vp. Vortices can move faster than fluxoids 
if LOi = LUcr, or they can move with the same velocity 
if —LOcr < oJi < ujcr, Or they move slower than fluxoids 
if LOi — —uJcr- To perform the numerical calculations, 
we start with Vp{Rc,tQ) — u„. This is a preliminary 
assumption, and we use it to calculate F^, F^, Fcmst, 
and Fn to find ui from equation (|40)) . Then, to find the 
actual velocity of the fluxoids, we check the following 
conditions. If —LOcr < i^i < i^cr, fluxoids and vortices 
comove, and we set Vp = w„. If w; > ojcr, the result 
derived from the preliminary assumption is not correct. 
Since, uoi can not be greater than lOc- Therefore, we set 
LOi — LOcn and calculate i^„(tt'cr)j where Vp is given by 
equation (|47)) . In this case vortices move faster than 



fluxoids. If uii < 



we set uji 



and calculate 



F„(— Wcr ). Equation (|47p gives us Vp, and therefore 7 
(see equation I53p . In this case vortices move slower than 
fluxoids. Having Vp , we calculate S{Rc,t) from equa- 
tion ()54p . This gives us the inner boundary condition 
for equation ([50)) . In our numerical calculations we use 
a normalized stream function s{r,t) = S(r,t) j {B^oRj^), 
where Bso is the initial value of surface magnetic field. 
The normalized stream function s has no units. Surface 
magnetic field in terms of the normalized stream function 
s is given by Bs{t) — Bsos{Rn,t). Therefore, from 
the core magnetic field is given by 



B,{t) = 2Bsos{r,t) 



(55) 



(56) 



Therefore, Fc is always negative. 



We assume that the initial profile of the normalized 
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TABLE II: Some trial values of the initial magnetic field of 
the neutron star, Bsi, intitial rotational (spin) period of the 
neutron star, Pai, the initial orbital period of the binary, Poi, 
the velocity of stellar wind, the rate of mass loss from a 
low main sequence companion, Af2, impurity parameter, Q, 
and the efficienty factor, is shown in the Table. 



(G) 


14, 
(kms~^) 


p 

(s) 


Log{M2) 
(Me yr-^) 


Q 


Ma 
(Me) 


c 




800 


1 


-13 


1 


1 


10 


3 X 10^^ 


600 


0.4 


-14 


0.1 


0.9 


1 




500 


0.1 


-15 


0.01 


0.8 


0.6 


3 X lO" 


400 


0.01 


-16 






0.1 


10^^ 












0.02 



(iii) Accretion phase: Accretion phase starts when 
Rco > Rm- In this phase some portion of the stellar wind 
is absorbed by the neutron star, therefore a is equal to 



1 



M2' 



(60) 



where Mn — —Ma- Finally, /? is given by equation 
Accreted matter, when it reaches the magnetosphere, in- 
teracts with the neutron star. In this phase Vd < 0, and 
the neutron star spins up according to equation ^2]) . The 
evolution of the magnetic field in the crust is governed 
by ([50)1 . where v is found from (|77|) . 



NUMERICAL RESULTS AND CONCLUSION 



stream function is given by 



s(r, to) = 1 - Sc 



Rn-r 

Rn — Rc 



(57) 



where ^ Sc ^ 1. Sc — 1 gives an initial profile with no 
magnetic field in the core, while Sc — gives a profile of 
uniform magnetic field. In terms of this initial profile the 
initial ratio of core to the surface magnetic field is given 
by 



— =2(1-S,) — 



(58) 



It is possible to consider that initially Bso = Bco- How- 
ever, this is not the only possibility, and early thermo- 
dynamic instabilities of the crustal layers might lead to 
Bso > BcO- As it was mentioned before we consider the 
evolution of a neutron star in different phases, 
(i) Ejector phase: In this phase the orbital separation 
changes according to the equation (fT8|) . where all the 
matter which is lost from the companion leaves the bi- 
nary, therefore a = f3 = 1. Since the stellar wind of the 
companion leaves the binary, the mass of the neutron star 
remains constant, = 0. Due to the magnetodipole 
radiation the neutron star spins down at the following 
rate 



dPs 
dt 



8 f BlRjsm'x 
c^I 



1 



(59) 



where x the angle between the rotational and magnetic 
axis. We set this angle to 7r/2. The crustal magnetic 
field evolution is governed by equation (j30p . where the 
material movement is set to zero, v = Q. Equation (j51|) 
gives the evolution of the stream function. 

{11) Propeller phase: In this phase we have the same set 
of equations, but instead of equation ([55]) we use ((T^)- 
> in this phase and the neutron star spins down. In 
this case a — 1 and (3 is given by equation P^ . Change 
in the mass of the neutron star is zero. In equation pop . 
the velocity of the material movement in the crust is zero. 
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FIG. 1: The evolution of the surface magnetic field for: (1) 
A binary neutron star with a low main sequence compan- 
ion, the effect of the superconducting core included. Initially 
the neutron star is in ejector phase. The first circle corre- 
sponds to the point of the transition to the propeller phase 
and the second circle corresponds to the point of the transition 
to the accretion phase. (2) A binary neutron star with the 
same companion, the effect of the superconducting core not 
included. Initially the neutron star is in ejector phase. The 
first cross corresponds to the point of the transition to the 
propeller phase and the second cross corresponds to the point 
of the transition to the accretion phase. (3)A single neutron 
star, the superconducting core is included. The binary pa- 
rameters are M2 =Mq, R2 =Rq, Poi = 14 days, M2 = 10"^^ 
Mq/jt, ( = 0.1, Vm = 500 km/s. Neutron star parameters 
are P,, = 0.01 s, B,, = 10^^ G, B„ = 10^° G, and Q = 0.01. 
Diamonds are observed pulsars with low mass main sequence 
companions {M2 = 0.8— 1.5 Me). At the given age the value 
of the surface magnetic field of these pulsars in known. 



We present the results of our numerical calculations for 
a neutron star of mass Af„ = l.SSM©, radius i?„ = 7.921 
km, and the radius of the core Rc = 7.43 km, where 
we assume that crust extends up to the density P c = 
2.4 X 10^** g cm"-^. We use the BPS equation of state [3l| 
in the crust of the neutron star. One need to test the 
results for various equations of states and compare them 
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FIG. 2: The evolution of the core magnetic field for: (1) the 
binary neutron star with a low main sequence companion, (2) 
the single neutron star. Parameters are the same as in Fig 1. 




FIG. 3: The evolution of the magnitude of the velocity of the 
neutron vortices \Vn\ and proton fluxoids |Vp|. The circle 
and the triangle mark the points when the velocities change 
sign and become negative. Parameters are the same as in Fig 
1. 



FIG. 4: The evolution of the magnitude of the drag force |Fi,| 
and pinning force |F„|. Circles mark the points when the 
forces change sign. F„ is negative after the circle and F„ is 
negative before the circle. Parameters are the same as in Fig 
1. 
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FIG. 5: The evolution of the magnitude of the buoyancy force 
I Ft I and |Fc|. The buoyancy force Fj, is always positive and 
Fc is always negative. Parameters are the same as in Fig 1. 



to the population of the MSPs. We use the cooling data 
from Neutron Star Theory Group at UN AM 34]. We do 
not take into account the dependence of the temperature 
into the densities inside the neutron star. For the further 
improvement of the model this effect has to be taken into 
account. Data of the observed pulsars is taken from [s^ 
(refer to Table III) . We assume a very small value for the 
moment of the inertia of neutron star / — 10^° g cm^^. 
Such of the moment of inertia is possible if the center of 
the neutron star is very dense and compact (36j (we have 
used BPS equation of state only as an approximation in 
the crust). This small value of the moment of inertia 
causes the fast rate of change of the spin period of the 
neutron stars with low mass main sequence companions 
in the propeller and accretion phases. We show that this 
scenario is in favor of the observations. Similar values of 
the moment of inertia have been assumed by other au- 
thors [5|]. However, the main features of this model are 



independent of the value of the moment of inertia. Nu- 
merical calculations were implemented in C-\ — h code. We 
consider some reasonable values for the initial magnetic 
field of the neutron star, B^i, intitial rotational period 
Psi, the initial orbital period of the binary Poi, speed of 
the stellar wind , rate of mass loss from the low main 
sequence companion M2, impurity parameter Q, and the 
efficiency factor C (see TABLE HH. We assume that M2, 
C, Vw (up to the Roche Lobe phase), and Q are constant 
parameters of the evolution. For simplicity, we assume 
that an accretion rate M2 = 1O"1*-1O"^°M0 yr"^ heats 
up the neutron star up to Tac = 3 x 10^ K, and an ac- 
cretion rate M2 = 10"^" — 10^^ M0yr~^ or higher heats 
up the neutron star up to Tac = 3 x 10^ K. However, the 
neutron star might enter the accretion phase at a rela- 
tively young age, when its surface temperature is higher 
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FIG. 6: The evolution of the surface magnetic field for the 
binary neutron star corresponding to different values of initial 
core magnetic field Bd- Parameters are the same as in Fig 1. 
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FIG. 7: The evolution of the surface magnetic field for a bi- 
nary neutron star corresponding to different values of initial 
surface magnetic fields Bs^. (1) B,i = 10^'' G. (2) B^, = 10^ 
G. Initial value of the core mag netic field is Bd = lO'' G. The 
rest of the parameters are the same as in Fig 1. 



than Tac- In this case, for simphcity, we let the neutron 
star cool down freely to the accretion temperature Tac, 
and then keep the star at this temperature as long as the 
accretion phase continues. 

We considered the evolution of a binary neutron star. 
Our aim was to find a more realistic model explaining 
the formation of MSPs, taking into account the effects of 
the core magnetic field and orbital evolution. The influx 
of the magnetic field is a unique signature of a binary 
neutron star with a superconductive core. In this model 
the magnetic field and the spin evolution play the most 
important roles. We were able to track the evolution of a 
neutron star in the Bg — Ps diagram. The spin, magnetic 
field, and the orbital evolution are coupled together. We 
consider only the low mass main sequence companions. 
The evolution has been considered for a period of 5 x 10^ 
yrs. The evolution has been stopped before the Roche- 



FIG. 8: The evolutionary track of surface magnetic field Bs 
and spin period for t = 4 x 10^ yrs. (1) Ah = 10"" 
M0/yr. The square shows the end state of the evolution after 
i = 5 X 10^ yrs. The triangle shows the end state of the 
evolution after t = 7 x yrs. (2)M2 = 10"^^ Mg/yr, The 
square shows the end state of the evolution after t = 5 x 10^ 
yrs. Parameters are: Psi = 0.01 s, Bsi = 5 x 10^^ G, Bd = 
10^^ G, Poi = 5 days, Q = 0.3, and ^ — 0.1. Black circles are 
observed millisecond pulsars. 
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FIG. 9: The evolutionary path of a binary neutron star in 
B — Ps diagram corresponding to different values of initial 
surface magnetic fields Bsi. (1) Est = lO" G. (2) Bsi = 10^° 
G. Initial value of the core magnetic field is Bd = 10^° G. 
The rest of the parameters are the same as in Fig 1. 



Lobe overflow. 

Fig [1] and [5] compare the evolution of a binary and 
a single neutron star, and the effect of the supercon- 
ductive core. For the single neutron star the magnetic 
field always decreases. For the binary neutron star if 
the effect of the core magnetic field is not included the 
surface magnetic field always decreases. For the binary 
neutron star with the superconductive core the magnetic 
field starts to increase at some point during the accretion 
phase. During the accretion phase the neutron vortices 
move inward to the core (Fig [3]). Inward motion of the 
neutron vortices does not imply the inward motion of the 
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FIG. 10: The evolution of the spin period of a binary neutron 
star. (1) Bsr = 5 X lO" G and M2 = 10"^^ Mo/yr. (2) 
Bs^ = 5x 10" G and Ah = 10"" Mo/yr. (3) Bs^ = xlO" 
G and M2 = 10"^^ Mo/yr. The initial value of the core 
magnetic field is Bd = 10^0 G and the initial orbital period 
is Poi = 4 days. The rest of the parameters are the same as in 
Fig 1. Triangles are observed binary pulsars. These pulsars 
have an orbital period 0.1 — 3.5 days at their current age. 
Their companions are low mass M2 = 0.8 — 1.3 Mq. 




FIG. 12: The evolution of the surface magnetic field for a 
binary neutron star corresponding to different values of M2. 
(1) M2 = lO-i^Mo/yr. (2) Ma = IQ-^^Me/yr. (3) M2 = 
lO~^^M0/t/r. Initial value of the core magnetic field is Bd = 
10^^ G, the initial value of the surface magnetic field is Bsi = 
10^^ G, the initial orbital period is Poi = 10 days, and Q = 
0.01. The rest of the parameters are the same as in Fig 1. 
Circles mark the point of the evolution where Ps = Peg for 
the first time. Magnetic field increases by a tiny amount after 
this point. 
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FIG. 11: The evolution of the surface magnetic field for a 
single neutron star corresponding to different values of initial 
core magnetic field Bd- Parameters are the same as in Fig 1. 



proton fluxoids. To derive the direction of the motion of 
the fluxoids we have considered various forces which act 
on the proton fluxoids at the core-crust boundary (Fig [4] 
and Fig [5]). Finally at some point of the evolution flux- 
oids start to move inward into the core (Fig [3]) . When 
the fluxoids start moving inward into the core, the core 
magnetic field increases. Be gives us the inner boundary 
condition for equation (|30p . The increase in the strength 
of the core magnetic field causes the enhancement of the 
surface magnetic field of the neutron star. One might 
wonder if the magnetic field can increase to such a large 
values that would destroy the standard scenario of the 
formation of the MSPs. Fortunately, this is not the case. 



A neutron star does not spin up to arbitrary short spin 
periods. There is a minimum limit on the spin period of 
a neutron star 1141 When this limit reaches the neutron 
star spin-up phase stops (Ps = 0) temporary as long as 
Ps = Peq (as long as the neutron star is on the spin-up 
line). Consequently Vn = as long as Pg — Peq- We 
assume that Vp = Q when w„ = 0. Therefore, the increase 
of the core magnetic field stops at this point, at least 
temporarily. If the core magnetic field is kept constant 
for some period of the time, the surface magnetic field 
starts decreasing again (refer to the short tail of the Fig 
[S]). To have another period of the spin-up, the surface 
magnetic field of the neutron star should decrease. This 
will cause the decrease of Peq- Therefore, the neutron 
star may enter another short period of the spin up and 
increase of the magnetic field. In this case the neutron 
star might spend a very large amount of time without 
moving much in the Bs — Ps diagram. Pulsars with a 
weak binary, low rate of mass loss from the companion 
M2 , need much more time to enter to the region of the 
MSPs. Fig [8] presents the evolution of a neutron star in 
the Bg — Ps diagram. Pulsar (1) in this figure with a com- 
panion that is losing mass at rate M2 = lO~"'^''M0yr~^ is 
very far from the region of the millisecond pulsars after 
t — A y. lO^yrs of the evolution. After 7 x lO^yrs of the 
evolution this pulsar joins the region of the millisecond 
pulsars. For a binary neutron star the final value of the 
surface magnetic field and the position of the neutron star 
in the Bg — Pg diagram is not very sensitive to the initial 
strength of the core magnetic field or the initial strength 
of the surface magnetic field (Fig . This is due to 
the existence of the minimum spin period for a neutron 
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FIG. 13: The dependence of the final surface magnetic field 
Bsf on the rate of the mass loss from the companion M2- 
ii)Poi = 4 days (2)Poi = 10 days. The initial strength of the 
surface magnetic fields is Bsi = 10^'' G, the initial strength of 
the core magnetic field is Bd = xlO^ G, Q = 0.01, C = 0.1, 
and the rest of the parameters are the same as in Fig 1. The 
evolution has been considered for 5 x 10^ yrs. Dash line marks 
the range of the strength of the surface magnetic field which 
is the most reasonable for millisecond pulsars Bs = 10^ — 10^° 
G. The lowest value of the surface magnetic field observed is 
6.67 X lO'^G which belongs to ,72229 + 2643. This pulsar is 
3.2 X 10^° yrs old. 
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FIG. 14: The dependence of the final surface magnetic field 
Bsf on the initial orbital period of the binary Poi for two 
different values of the speed of the stellar wind: (1)14; = 500 
km/s {2)Vw = 800 km/s. The initial strength of the surface 
magnetic fields is Bsi = 10^^ G, the initial strength of the core 
magnetic field is B„ = xlO^^ G, Q = 0.01, C = 0.1, and the 
rest of the parameters are the same as in Fig 1. The evolution 
has been considered for 5 x 10^ yrs. Triangles are observed 
pulsars. These pulsar have an age 8.78 x 10* - 4.83 x 10^ 
yrs, and low mass companions M2 — 0.8 — IMq. Their initial 
orbital period have been approximated based on their current 
orbital period within our code. Squares mark the point of 
the transition from the ejector to the propeller and from the 
propeller to the accretion phase. 



FIG. 15: The evolution of the surface magnetic field for: 
(1)Q=0.3 (2) Q=0.001 (3)Q=1. The initial surface magnetic 
filed is Bsi = 5 X 10^"^ G and the initial core magnetic field is 
Bci — 10^^ G. Parameters are the same as in Fig 1. 
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FIG. 16: The evolution of the surface magnetic field for a 
binary neutron star. (1) The orbital evolution is taken into 
account. (2) The orbital evolution is neglected. Initially the 
neutron star is in ejector phase. The first cross corresponds 
to the point of the transition to the propeller phase and the 
second cross corresponds to the transition to the accretion 
phase. Presence of the orbital evolution does not effect the 
transition points much. Parameters are: M2 = 10"" Mo/yr, 
Bsi = 10^^ G, Bci = 10" G, Poi = 4 days, and the rest of the 
parameters are the same as in Fig 1. 



star Peq. Neutron stars with different initial strength of 
the core or surface magnetic field reach Peq at almost the 
same value of the surface magnetic field (the same value 
of Peg), ahhough at different times. Fig [TUl shows that 
two pulsars with different initial surface magnetic fields 
reach the same value of the spin period. After this point 
the surface magnetic field remains almost a constant for 
10^ — 10^ yrs before it starts decreasing again. This fea- 
ture is unique to binary neutron stars. For example, the 
final stage of single pulsars is very sensitive to the ini- 
tial strength of the core magnetic field (Fig [TT|l . On the 
contrary, the value of the final surface magnetic field is 
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FIG. 17: The evolution of the orbital period Po for Poi = 4 
days and Poi = 4.6 days. B^i = 5 x lO" G, B„ = lO" G, 
M2 = lO'^^Mo/yr, and the rest of the parameters are the 
same as in Fig 1. 



very sensitive to the parameters of the binary, such as 
the rate of the stellar mass loss from the companion M2 , 
the initial orbital period of the binary Poi, the speed of 
the stellar wind V^, and the coefficient factor (. Neu- 
tron stars with different values of M2 reach the minimal 
spin period Peg at different values of the surface mag- 
netic field (Fig fT2|) . Therefore, they end up at different 
values of Bsf ■ Fig [5] shows that also the path of a neu- 
tron star in B — Ps diagram strongly depends on M2- 
Fig [13] illustrates the dependence of the finial strength 
of the surface magnetic field to the rate of the mass loss 
from the companion AI2 for two different values of the 
initial orbital period. Note that not all these pulsars 
have reached their Pgq in the period of the evolution we 
have considered, i = 5 x lO^yrs. The dashed line marks 
the most reasonable range of the surface magnetic field 
of the observed MS pulsars Bs = 10^ - lO^^G. Some of 
the pulsars which are below this line will gain a stronger 
surface magnetic field at later stage of their evolution. 
This feature has been illustrated in Fig[9l This property 
of these pulsars is due to the low value of M2. Namely, 
these pulsars need more time to reach their minimal spin 
period Peg- Fig [Til illustrates the dependence of the final 
strength of the surface magnetic field of a neutron star 
on the initial orbital period of the binary for two differ- 
ent values of the speed of the stellar wind Therefore, 
this figure also shows that the evolution depends on Vw 
Pulsars in binaries with long orbital period need to spend 
more time in the accretion phase to end up in the region 
of the millisecond pulsars. The magnetic field increases 
during this phase and the spin period decreases. Even- 
tually, the spin period reaches its minimal value at high 
enough magnetic field. At this stage the pulsar reaches 
the region of the millisecond pulsars. 

The evolution of the surface magnetic field is very sen- 
sitive to the different values of Q, due to the dependence 
of the impurity conductivity on this parameter fFig llSp . 



FIG. 18: Initial and final state of binary neutron stars in 
B ~ Pa diagram for difi^erent trial parameters. 5 x 10^ yrs of 
the evolution has been considered. Gray circles are observed 
pulsars. MSPs are located at the lower left region of the dia- 
gram (region of the low magnetic field and short spin period). 
Not all the neutron stars with different parameters end in the 
MSPs region in 5 x 10® yrs. (1) Initial point of the evolution of 
a young neutron star. Bsi — 10^* G, Psi ~ 0.01 s. (la) Final 
state of the evolution of the neutron star (1) for M2 — 10~^^ 
Mo/yr, C ~ 0.1, Poi = 4 days, (lb) Final state of the evolu- 
tion of the neutron star (1) for M2 = 10"" Mo/yr, C = 0.1, 
and Poi = 4 days. (Ic) Final state of the evolution of the 
neutron star (1) for M2 = 10"" Mo/yr, C = 0.1, and Poi = 4 
days. (Id) Final state of the evolution of the neutron star (1) 
for M2 = 10"^^ Mq/ yr, C = 0.1, and Poi = 4 days. (If) Final 
state of the evolution of the neutron star (1) for M2 ~ 10"^^ 
Mo/yr, C = 0.1, and Poi = 100 days. (2) Initial point of the 
evolution of a young neutron star. Bsi = lO" G, Psi = 0.004 
s. (2a) Final state of the evolution of the neutron star (2) 
for A/2 = 10"^^ M©/yr, ( = 0.1, and Poi = 4 days. If one 
considers the evolution of (1) or (2) of a single neutron star 
the end state is Bsf — 10® G and Psf = 2.34 s. Pulsars (2a), 
(Ic), (If), and (Id) are still in the graveyard at t — 5 x 10® 
yrs. 



We have not considered different equations of state and 
we have assumed a uniform temperature inside the neu- 
tron star. However, the dependence of the surface mag- 
netic field on Q implies a sensitive dependence of the 
evolution on the conductivity. The conductivity depends 
on the density at different layers inside the neutron star 
(equation of state) and on the temperature. Therefore, 
various equations of the state have to be tested. 

FigfTHlillustrates that the effect of the orbital evolution 
in very close binaries can be very important at the very 
late stages of the evolution. Fig [17] shows the evolution 
of the orbital period of the binary. 

In this model the magnetic field and spin period of the 
pulsar change substantially. Therefore, binary pulsars 
spend some large period of their life in the graveyard (Fig 
m [9l and [T8)) . Pulsars with weaker binaries remain for 
longer periods in the graveyard. Perhaps this accounts 
for the very low population of the MSPs in the regions 
belonging to the normal pulsars. We propose that in 
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FIG. 19: The evolutionary track of the surface magnetic field 
Bs and spin period Ps for t = 4x 10^ yrs. Two different values 
of the moment of the inertia of the neutron star is chosen. (1) 
/ = 10*° g cm^ (2) / = 0.3 X 10"*^ g cm^ The rest of the 
parameters are the same as Fig 1. 



the moment of inertia. 
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TABLE III: Example of the observed pulsars with low mass 
companions. Most of these companions are neutron stars 
(NSs) or white dwarfs (WDs). Unfortunately, there are not 
many observed pulsars with low mass main sequence (MS) 
companions. This is not the whole set of observation data 
we have used. This data for this table is taken from [35{ . 
We have rounded the numbers here (in figures more accurate 
numbers are used). 

Name CS Ps Po Ma Age 

(s) (days) (yrs) (G) 



other models of the evolution where the effect of the su- 
perconductive core is neglected could not be as successful 
as the current model in describing the absence of the bi- 
nary pulsars from the region of normal pulsars. In these 
models ([1], magnetic filed always decreases and final 
value of the magnetic field should be between lO'^G and 
10^"G. If pulsars in such models get to very low values 
of the surface magnetic field they can not come back to 
the region of MSPs. 

In Fig [19] we compare the evolution of two pulsars of 
small and large moments of inertia. Pulsar (2), with large 
value of the moment of inertia, enters the propeller phase 
at age 3 x 10^ yrs, and later enters the accretion phase 
at age 9 x 10^ yrs. However, due to the very large value 
of the moment of inertia its spin does not change much. 
Therefore, this pulsar follows almost vertical paths in the 
Bs — Ps diagram. However, this pulsar should be visible 
at the propeller and accretion phases. If this were true, 
we would observe many binary pulsars at different ages, 
passing through the region of normal pulsars. Depending 
on the different values of their spin period and magnetic 
field, we would find them at different paths in Bs — Ps 
diagram. Almost all of these pulsars would pass almost 
vertically through the region of normal pulsars. Accord- 
ing to the observations, there are only a handful of binary 
pulsars in the region of normal pulsars (only 2%). There- 
fore, one needs an evolution that explains the absence of 
binary pulsars from the region of normal pulsars. This 
type of the evolution can be provided by low values of 



0751 + 1807 


MS 0.003 


0.26 


0.1 


7.0 


X 


10^ 


1.67 X 10* 


J1903 + 0327 


MS 0.002 


95.17 


0.89 


1.8 


X 


10^ 


2.0 X 10* 


J1913-I-16 


MS 0.05 


0.32 


0.8 


1.0 


X 


10* 


2.28 X 10^° 


B1718 - 19 


MS 1.0 


0.2 


0.1 


9.8 


X 


10** 


1.29 X 10^^ 


J0621 4- 1002 


WD 0.03 


8.3 


0.4 


9.6 


X 


10^ 


1.1 X 10^ 


J2145 + 0750 


WD 0.01 


6.8 


0.4 


8.5 


X 


10^ 


7.0 X 10* 


J1435 - 6100 


WD 0.01 


1.3 


0.9 


6.0 


X 


10^ 


4.8 X 10* 


J1157 - 5112 


WD 0.04 


3.5 


1.18 


4.8 


X 


10^ 


2.5 X 10^ 


J1528 - 3146 


WD 0.06 


3.2 


0.9 


3.9 


X 


10^ 


3.9 X 10^ 


J1802 - 2124 


WD 0.01 


0.7 


0.8 


2.8 


X 


10^ 


9.7 X 10* 


J1454 - 5846 


WD 0.05 


12.4 


0.86 


8.8 


X 


10* 


6.2 X 10® 


J0737 - 3039.4 


WD 0.02 


0.10 


1.3 


2.0 


X 


10* 


6.4 X 10® 


B2303 + 46 


WD 1.06 


12.3 


1.1 


2.9 


X 


10^ 


7.8 X 10" 


J1141 - 6545 


WD 0.39 


0.2 


0.9 


1.5 


X 


10" 


1.3 X 10^^ 


B1820 - 11 


WD 0.39 


0.2 


1.0 


1.4 


X 


lO** 


1.3 X 10^^ 


J1022 + 1001 


0.01 


7.80 


0.7 


6.0 


X 


10^ 


8.5 X 10* 


B0655 + 64 


0.19 


1.02 


0.6 


4.5 


X 


10^ 


1.1 X 10^" 


J1439 - 5501 


0.03 


2.12 


1.11 


3.2 


X 




2.0 X 10® 


B1820 - 11 


0.09 


13.6 


0.4 


1.5 


X 




9.7 X 10® 


J1753 - 2240 


0.27 


357 


0.6 


3.2 


X 


10" 


6.2 X 10" 


J1748 - 2446 


0.008 


0.38 


0.5 










J1811 - 1736 


NS 0.10 


18.78 


0.9 


1.8 


X 


10"' 


3.1 X 10® 


J1829 + 2456 


NS 0.04 


1.18 


1.26 


1.2 X 10"' 


1.5 X 10® 


J1756 - 2251 


NS 0.03 


0.32 


1.0 


4.4 X 10* 


5.4 X 10® 


B1534 + 12 


NS 0.04 


0.42 


1.3 


2.5 X 10* 


9.7 X 10® 


B1913 4-16 


NS 0.06 


0.32 


0.8 


1.1 


X 


10* 


2.3 X 10^" 


J1906 -f 0746 


NS 0.14 


0.16 


0.8 


1.1 X lO'^ 


1.7 X 10^^ 



Natural Science and Engineering Research Council of 
Canada. 



[1] Srinivasan G., Bhattacharya D., Muslimov A. G. and (2001). 

Tsygan A. J., Current Science , 59, p. 31 (1990). [4] M. Jahan-Miri, ApJ, 532, p. 512 (2000). 

[2] Konenkov D., and Geppert U., A&A, 372, p. 583 (2001). [5] Urpin V., Geppert U., and Konenkov D., MNRAS, 295, 
[3] Konenkov D., and Geppert U., MNRAS, 325, p. 426 p. 908 (1998). 



16 



[6] Ding K. Y.,Cheng K. S. , and Chau H.F, ApJ,408, p.l67 

(1993). 

[7] Jahan Miri M. and Bhattacharya D., MNRAS, 269, p. 
455 (1994). 

[8] Pringle J.E and Rees M.J, A&A, 21, p.l (1972). 
[9] Illarionov A.F and Syunyayev R.A, A&A, 39, p. 185 
(1975). 

[10] Thompson C. and Duncan R.C, ApJ, 408, p.l94 (1993). 

[11] Jahan Miri M, ApJ, 532, p. 514 (2000). 

[12] Van Riper K. A., ApJS, 75, 449 (1991). 

[13] Hunt R., MNRAS, 151, p. 141 (1971). 

[14] Wang Y. M., A&A, 102, p. 36 (1981). 

[15] Haensel R, Zdunik J. L., A&A, 227, 431 (1990). 

[16] Fujimoto M., Hanava T., Icko Iben Jr., Richardson M., 
ApJ, 287, 813 (1984). 

[17] Note that it is still possible to have a neutron star in the 
ejector, or the propeller phase when the Roch lobe over- 
flow starts. However, this happens only if a young neu- 
tron star has an old companion. This is because young 
neutron stars have a strong radiation pressure. These 
cases are rare and we do not consider them. 

[18] Paczynski B., "Evolutionary processes in close binary 
systems" (1971). 

[19] Hoyle F. and Lyttleton R.A, POPS, 35, p. 405 (1939). 

[20] Abbott D., ApJ, 225, p. 893 (1978). 

[21] Ostriker J.P and Davidson K., lAUS, 55, p. 143 (1973). 

[22] Van den Heuvel E. P. J., Ann. NY. Acad. Sci, 302, p.l4 



(1977). 

[23] Verbunt F. and Zwaan C, A&A, 100, p. L7 (1981). 

[24] Smith M. A. , PASP, 91, p. 737 (1979). 

[25] Abdolrahimi S. " Dynamics of Angular Momentum 

Transfer into neutron stars from Disk Accretion", Ms. 

Thesis, Alzahra University (2006). 
[26] Bhattacharya A. D. and Datta B., MNRAS, 282, p. 1059 

(1996). 

[27] Konar A. S. and Bhattacharya A. D., MNRAS, 303, 
p.588 (1999). 

[28] Geppert U. and Urpin V., MNRAS, 271. p. 490G (1994). 
[29] Urpin V.A, and Yakovlev D.G, SvA, 24, p. 425 (1980). 

[30] Muslimov A. and Tsygan A., Ap&SS, 115, p. 41 (1985). 
[31] Goldreich and Rcisencggcr, ApJ, 395, p. 250 (1992). 
[32] Jones P. B. , MNRAS, 233, p. 875J (1988). 
[33] Baym G., Pcthick C., and Sutherland P., ApJ, 170, 
p.299(1971). 

[34] Buccheri R., Paradijs J.V. , and Alpar M. A. , Obs, 119, 

no.ll49, p.105 (1999). 
[35] ATFN pulsar catalogue. Manchester R. N., Hobbs G. 

B., Teoh A., and Hobbs M., AJ, 129, 1993 (2005), 

(www.atnf.csiro.au). 
[36] J. M. Lattimer and M. Prakash, ApJ, 550, p. 426 (2001). 
[37] Nagase F., PASJ, bf 41, p.l (1989). 
[38] Illarionov A. F. and Sunyaev R., A&A, 39, p. 185 (1975). 



